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Jet production in electron-proton collisions at HERA provides a unique testing 
ground for Quantum Chromodynamics (QCD). Apart from the determination of 
the strong coupling constant Os, ep jet data may especially be used to gain insight 
into the dynamics of the exchanged parton cascade, whose structure is probed by 
the high-iJ-ji dijet system; thus information on the parton content of the proton 
and (quasi-)real and virtual photons is obtained. This report touches some of these 
aspects revealed in recent jet data from the HERA experiments which are testing 
perturbative QCD at the limits of applicability. 

1 Introduction 

Multi-jet production in deep inelastic ep-scattering (DIS) provides special 
sensitivity to the mechanisms of the strong interaction and can hence be used 
to test the predictions of perturbative Quantum Chromodynamics (pQCD) in 
a rather inimitable way. This is especially true since the HERA experiments 
deliver data over a large range of both the four-momentum transfer, Q^, of 
the exchanged photon and the transverse energy of the observed jets, Et- 

In lowest order, multi-jet production in DIS is described by the QCD- 
Compton and the boson-gluon fusion processes (Figure 0), with the momen- 
tum distributions of the incoming partons taking part in the hard interaction 
given by the parton density functions (PDFs) of the proton; the Q^-evolution 
of the latter is described by the DGLAP equationsEJ which in lowest order 
are equivalent to the assumption of exchanging a strongly fct-ordered parton 
cascade, fcf being the transverse momentum of the partons within the cascade. 
However, as illustrated by the generic diagram in Figure ||c, to calculate ep 
jet cross sections requires additional terms, such that in certain parts of phase 
space, where fej-ordered parton emission is no longer manifest, the standard 
DGLAP approach has to be extended e.g. with the concept of photon struc- 
ture, which preserves the perturbative ansatz of the DGLAP evolution scheme 
but results in an artificial violation of the fct-ordering. Within this concept, 
high Et jets can be produced not only by direct processes in which the vir- 
tual photon interacts as a point-like particle with a parton out of the proton, 
but also by resolved processes where the photon interacts hadronically. In 
the photoproduction limit (Q^ 0), the cross section is thus sensitive to the 
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Figure 1. Schematic Feynman diagrams for tiie dijet production meclianism witii (a) the 
QCD-Compton and (b) the boson-gluon fusion process; (c) generic picture for dijet produc- 
tion via photon-parton fusion where the actual interaction mechanism is unknown and to 
be resolved by the dijet system. The two scales involved in these processes are indicated as 
the four-momentum transfer and the transverse jet energy Et- 

parton distributions of both the photon and the proton at a scale set by the 
transverse energy of the jets, Et- 



The two scales involved in ep jet production, and Et, can now be 

used to resolve the structure of the photon-parton interaction in general. Sev- 
eral kinematic regions can be distinguished, depending on the absolute and 
relative sizes of these two scales. In the regime where and E^ are both 
large, perturbative methods are clearly justified and the DGLAP ansatz can 
be used to extract information on the proton structure and the strong coupling 
as- Typically, leading-order (LO) Monte Carlo models approximate higher or- 
der pQCD contributions in this regime by parton showers. If ^ E^ the 
jet cross sections become sensitive to the hadronic structure of the photon. 
The region where « E^ is of special interest, since in this regime the 
DGLAP ansatz is expected to break down for small values of the variable xb 
(Bjorken-x) revealing effects of non-fct-ordered parton cascades. Characteris- 
tics expected from other evolution schemes like BFKLH or CCFAfi could be 
observable in this phase space regime. 

Any deviations from next-to-lcading order (NLO) predictions may, how- 
ever, simply indicate that fixed order calculations beyond NLO are needed. 
These may also help to resolve other ambiguities associated with high Et jet 
production in DIS: where and E^ both are substantially large ^qcd) 
there is ambiguity in the choice of the renormalization scale^ It is important 
to know where these ambiguities have a significant effect on the cross sections. 

2 Testing pQCD at Large Scales 

Jet production at large and large Et has been studied by both the HI 
and the ZEUS collaborations. In this regime the choice of the renormaliza- 

"Analogously such ambiguity exists also for the factorization scale. However, since the 
effects due to the choice of the renormalization scale are generally larger, this ambiguity is 
not explicitly mentioned in the context of this paper. 
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tion scale {fir = vQ^j -^t) is of minor importance since scale dependencies 
of the cross sections are generally small if /i^ is large. Hence NLO order cal- 
culations are expected to provide a good description of the data as can be 
seen from Figurc.|^ which shows the inclusive jet cross section (Pa/dErdQi. 
measured by HlH. Jets have been defined using the inclusive k± algorithrra 
in the Breit frame where the photon collides head-on with the incoming pro- 
ton. Over the whole analyzed phasg space the DISENT NLO calculation!, 
using the CTEQ5M1 PDFs as inputll and corrected for hadronization effects 
(< 10%), agrees perfectly well with the data. The size of the scale uncertain- 
ties together with the error contribution from hadronization corrections can 
be inferred from Figure || (right). From these data, HI extracts the strong 
coupling constant as in bins of Et, chosen to be the renormalization scale. 
The QCD predictions are fitted to the jet cross sections using the CTEQ5M1 
parameterization of the PDFs and the strong coupling constant as the single 
free parameter. The result is shown in Figure ^ and clearly indicates the run- 
ning of as according to the renormalization group equation. A combined fit 
to all 16 data points shown in Figure ^ results in 

as{Mz) = 0.1186 ± 0.0030 (exp.)t°:0039 (theo.)t°:°°i (pdf) (m^ = Et), 

where the largest contributions to the experimental error come from uncer- 
tainties of the hadronic energy scale. The theoretical error is dominated by 
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Figure 2. The left figure shows the inclu- 
sive jet cross section as a function of the 
transverse jet energy Et in different regions 
of Q^. Jets are reconstructed in the Breit 
frame using the inclusive kj_ algorithm. The 
data are compared to the NLO pQCD pre- 
diction with and without hadronization cor- 
rections. For two of the Q^-bins the right 
figure shows the corresponding ratio of data 
to the theoretical prediction; the uncertain- 
ties due to hadronization effects and the 
choice of renormalization and factorization 
scales are indicated by the band. 
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Figure 4. The correlation of as and the 
gluon density in the proton from a fit to 
HI jet and inclusive DIS data for diff'er- 
ent values of the fractional momentum 
X of the gluon. The central fit result is 
indicated by the full marker while the 
error ellipses show the experimental and 
theoretical uncertainties. 
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Figure 3. Determination of as from the inclu- 
sive jet cross section for the renormalization scale 
fir = Et- The results are shown for each Et 
value (circles) including experimental and theoret- 
ical uncertainties. The single values are extrapo- 
lated to the Z^-mass (triangles). The final re- 
sult as{Mz) (rightmost triangle) is obtained in a 
combined fit. The curves represent the combined 
result and its uncertainties once at the Z^-mass 
(lower curves) and once evolved according to the 
renormalization group equation (upper curves). 
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the uncertainty in the hadronization corrections and the dependence on the 
choice of the renormahzation scale. The uncertainty in the knowledge of the 
parton density functions has been estimatedu using the correlated errors pro- 
vided by a recent QCD analysisO together with its global PDF-fits. This error 
contribution is substantially reduced when using the dijet rate i?2-i-i to extract 
tts as done by the ZEUS coUaboratioro, leading to the preliminary result 



0.1166 ll]:E;EJf7(exp. 



+0.0055 
-0.0042 



(theo.);°:™i? (pdf) {^ir 



Both measurements of as depend on the knowledge of the parton con- 
tent of the proton. Vice versa, by using the best knowledge on one can 
use (multi-)jet cross section measurements to extract the quark and gluon 
densities in the protonQ. A test of pQCD independent of data from other 
experiments, has been made by the HI coUaborationQ in performing a simul- 
taneous determination of both quantities. This is done in a fit to the inclusive 
jet and dijet cross sections together with the inclusive DIS cross section, where 
the latter is restricted to the kinematic range 150 < < 1000 GeV^ and only 
constrains the quark densities in the proton. The results of this fit are shown 
in Figure^ as a correlation plot between as{Mz) and the gluon density xg{x) 
evaluated at four different values of the gluon momentum fraction x = 0.01, 
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0.02, 0.04 and 0.1. While the present data do not allow a simultaneous deter- 
mination of both parameters with competitive precision the sensitivity to the 
product as • xg{x) can clearly be seen. 

3 Scale Ambiguities 

The relevance of the choice of the renormalization scale at smaller can be 
seen from Figure ||, showing the dijeLxross section as a function of log as 
measured by the ZEUS coUaborationEj. As for the inclusive jet cross section, 
pQCD in NLO describes the data down to values of w 150 GeV^. However, 
at values below 150 GeV^ scale uncertainties become large and the choice 
of the renormalization scale is no longer irrelevant: while NLO calculations 
using /ir = v^Q^ describe the dijet cross section down to k. 10 GeV^, this 
is no longer the case for /i^ — Et- _ 

This issue has also been studied by the HI coUaborationO investigating di- 
jet event rates, i?2, at low values of the Bjorken-x variable 10^^ < xb < 10^^ 
and low Q^, 5 < < 100 GeV^. For different requirements on the trans- 
verse energies £'r,(i,2) of the final state jets, NLO QCD calculations are 
confronted with the data using two different choices of the renormalization 
scale = \/ and fi^ = ^/Q"^ + E^. Figure |^ summarizes the prelimi- 
nary result of this analysis in two representative {xB,Q^)-hms. At large 
{Q^ — 71 GeV^) and large xb (xb — 4.7 ■ 10^'^) the theoretical calculation is 

Figure 5. (a) The dijet cross section 
da / dlogQ^ for jets in the Breit frame 
using the inclusive fcx algorithm. The 
points represent the data with statisti- 
cal (inner bars) and total (outer error 
bars) uncertainties. The shaded band 
indicates the systematic uncertainty 
due to the energy scale of the ZEUS 
calorimeter. While the theoretical pre- 
dictions using /ir = \/ describe 
the measured cross section rather well, 
the dashed-dotted line representing a 
NLO calculation with /^r = 53 i/2, 
where Et^i are the transverse jet en- 
ergies of the jets, falls below the data 
for Q2 < GeV^. This can also be 
seen from (b) which shows the ratio of 
the measured cross section to the NLO 
prediction_using CTEQ4M PDFs and 
Mr = 
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Figure 6. Dijet rate R2 at two 
values of xg, as a function 
of the transverse energy require- 
ment Et^i > (5 + A) GeV on 
the jet with larger Et; for the 
second jet the transverse energy 
i?T 2 rnust be above 5 GeV. 



rather insensitive to the choice of the scale (Fig. ||, right), and — in agreement 
with the results from the dijet cross section measurement at high — gives 
a good description of the dataR In contrast, one obtains no safe theoretical 
prediction at low xb, (Fig.^, left): although data and NLO QCD calcu- 
lation do agree for fj,r = the large scale uncertainties lead to almost 
vanishing predictive power of the theory; choosing, however, a larger scale 
~ ^/Q"^ + E'^ such that scale uncertainties become small, leads to clear 
disagreement between the measured event rate R2 and the next-to-leading 
order QCD calculation, independent of the requirement on the transverse en- 
ergy of the jets Et,i > (5 + A) GeV. This indicates the necessity to include 
higher-order contributions by introducing, for example, the concept of photon 
structure. Choosing the same scale rtr-i= \/Q'^ + E^: comparison of the data 
with the NLO calculations by JetViPtj including a resolved-photon contribu- 
tion, indeed shows better, although not perfect agreementlld. 



4 Probing Parton Dynamics 

In order to further investigate the interplay of the two scales ^/Q^ and Ex-, 
the ZEUS collaboration has studiedo the cross section for jets produced at 
large pseudo rapidities as a function of the ratio E^/Q^. When compared to 
LO Monte Carlo models (Figure 0) the data are well described for E^ <^ Q^. 
However, onh^ those models which include non-fct-ordered parton emission 
(ARIADNE0, RAPGApH) reproduce the E^/Q'^ w 1 region; the fuU range 
can be described solely by models which include a resolved photon contri- 
bution as the RAPGAP Monte Carlo and a NLO calculation by JetViP. A 



''This is true unless a symmetric trapsverse energy cut is applied where the NLO QCD 
calculation becomes infrared sensitivd^^. In this region resummed calculations are needed 
but not yet available. Thus the presented measurement also provides an important reference 
for improved theoretical predictions. 
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Figure 7. Forward-jet cross section as a function of E^/Q^ in comparison (a) to several 
LO Monte Carlo predictions and (b) to RAPGAP and a NLO JetViP calculation both 
including resolved photon contributions. The shaded band corresponds to the energy scale 
uncertainty of the ZEUS calorimeter. 

reasonable description of these data is also found in a recent comparison0 
with Monte Carlo predictions based on the CCFMcl evolution equation which, 
by means of angular-ordered parton emission, is equivalent to the BFKLa 
approach for a; ^ while reproducing the DGLAP equations at large x. 

To explore a possible signature of BFKL dynamics, the w regime 
has been analyzed in more detail studying the Xs-dependence of forward par- 
ticle productiorj^. Figure || shows the a: ^dependence of the forward-jet and 
forward-7r° cross sections in this regionMlla in comparison to Monte Carlo 
models with and without a resolved photon contribution and — in case of 
the tt" cross sections — a LO BFKL calculation. The need for an additional 
contribution to the direct 7p-interaction can clearly be seen from these figures 
as the data are rather well described by models including a resolved photon 
component. Moreover, the ir^ measurement, which accesses very small poLax 
angles, is actually best described by a BFKL-based leading-order calculationl3 
when taking hadronization corrections into account. However, large theoreti- 
cal scale uncertainties diminish the significance of these comparisons. 



5 Jets in Photoproduction 

The production of hard dijet events in photoproduction is dominated by re- 
solved photon processes in which a parton in the photon with momentum 
fraction x-y is scattered from a parton in the proton. Hence studies of dijet 
production can be used to investigate the hadronic structure of the photon 



■^Particles produced at very small polar angels w.r.t. the incoming proton direction. 
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Figure 8. xs-dependence of forward- 
jet (left) and forward-Tr" (right) cross 
section measured for Ki Q^. The 
data are compared to Monte Carlo 
models with and without resolved pho- 
ton contribution as well as a LO BFKL 
calculation (right figure only). 



at « 0. In leading order the cross section is proportional to the pho- 
ton flux /^/g and maybe approximated by use of effective parton densities 
for the proton and photon such that adijet ~ f-y/e ' feB,-yfeS,p ■ where 
Ai is the matrix element of the hard parton-parton scattering process and 
/eff,[7,p] = /q/[7,p] + /q/[7,p] i^^/^ /g/[7,p]- With this relation HI extracts the 
gluon density of the photoncJ in leading order pQCD using the knowledge of 
/off,p from DIS and from the photon structure function F2 as mea- 

sured in e^e~ collisionsEll. The result is shown in Figure ^ indicating a strong 
rise towards low x~/. 

The determination of the gluon content of the photon is however strongly 
affected by non-pertutbative effects and the treatment of background from 
soft underlying eventgH3; it therefore can only be done using Monte Carlo 
programs which include phenomenological models to account for these effects 
but up to now are stilLrestricted to. leading order matrix elements. Other 
dijet analyses from HlE2l and ZEUScJ suppress these non-perturbative effects 
by harder Et requirements for the jets and thus investigate the parton dis- 
tributions in the photon at high x-y {x-y > 0.2) where quark contributions 
dominate]^. A comnaxison of the NLO prediction with the dijet cross section 
measured by ZEUSc3 for Er.jet > 14 GeV as a function of the pseudo ra- 
pidities of the two most energetic jets, 7]^"^* and 772°*, is given in Figure |lO| . 
While the NLO calculation describes the data at large x-y > 0.75 rather well, 



"^Actually the dependence on 



[ET,ie 



is studied, where Ee 



is the electron beam energy. In the leading order massless approximation x"^"" is equivalent 
to X-v. 
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Figure 9. The gluon distribution of tiie pho- 
ton multiplied by a~^Xj as a function of x^y 
at a mean scale of = 74 GeV'^. 
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a clear discrepancy is seen for the full a;^-range if both jets are going in the 
forward direction (positive 7]^°^). This is where substantial contributions from 
resolved processes {x^ < 0.75) are expected suggesting that, in the kinematic 
region of the measurement presented, the available paxameterizations of the 



6 Conclusion 

Measurements of ep jet production rates have been presented in a wide range 
of the four-momentum transfer and the transverse jet energy Et to study 
the dynamics of the underlying partonic interaction. At large scales fir = 



a/Q^, Et the data are successfully described by perturbative QCD allowing 

a determination of the strong coupling constant as and the gluon density in 
the proton. Scale ambiguities have been studied and, as gets smaller, 
sensitivity of the NLO predictions to the choice of the renormalization scale 
is observed. At low insight into the structure of the photon and effects of 
non-fcj-ordered parton evolution is revealed. 
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